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(54) Automatic control system which optimises eiectrophoretic separation 
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a knowledge base derived from calibration data obtained by separating macromoiecules, in particular DNA, of known size 
under a variety of planned conditions chosen to optimize separations. The user enters the smallest and largest 
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The system preferably uses the CHEF (Contour clamped Homogeneous Electric Field) system with a hexagonal 
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AUTOMATIC CONTROL SYSTEM FOR IMPROVED ELECT ROPHORESIS 
BACKGROUND OF THE INVEyrrXpfl 
This invention relates generally to gel 
electrophoresis and, more particularly, to a system for 
controlling contour-clasped homogeneous electric t\eld 
electrophoresis . 

The standard method used to separate, identify, 
and purify DNA and other macromolecules is 
electrophoresis through agarose gels. This technique is 
capable of efficiently and rapidly resolving mixtures of 
DNA fragments that cannot be separated adequately by 
other purification technologies. A gel containing a DNA 
sample is subjected to an electric field. This causes 
the individual macromolecules in the DNA to migrate 
through the gel. Because different DNA fragments migrate 
at different rates, visible "bands" of DNA form in the 
gel. The location of individual DNA fragments within the 
gel can be read directly from the bands. 

A particular DNA 1 s characteristics, including 
its size, charge, and conformation, influence its rate of 
migration and, hence, the separation that may be 
achieved. For example, molecules of linear and duplex 
DNA travel through gel matrices at rates that are 
typically inversely proportional to the log 10 of their 
molecular weights. Nevertheless, DNA macromolecules of 
identical molecular weight but different conformations, 
such as closed circular, nicked circular, and linear 
forms, migrate through agarose gel at different rates. 

In addition to DNA characteristics, 
electrophoretic parameters also influence the rate of 
migration (mobility) of macromolecules. DNA fragments of 
a given size migrate at different rates through gels 
containing different concentrations of agarose. By 
varying the concentration, it is possible to resolve a 
vide range of sizes of DNA fragments. The strength of 
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the applied field (voltage gradient), the ionic strength 
of the buffer, and the density of the superhelical twists 
m the closed circular form of DNA also influence 
migration. At low voltages, the mobility (cm/hr) of 
linear DNA fragments is proportional to the voltage 
applied. However, this linear relationship fails at 
higher potentials for high molecular weight fragments. 
Thus, the effective range of separation for a given 
agarose gels decreases as the voltage increases. 

Improved separation of macromolecules may be 
achieved by optimizing certain parameters. For example 
pulsed-field gel electrophoresis ( PFGE ) , a recent 
advancement, allows for greater resolution of DNA 
fragments for a given voltage. Schwartz and Cantor teach 
that lt is possible to resolve large DNA fragments in an 
agarose gel subjected alteratively to two alternately 
pulsing electrical fields at about 90-. (Schwartz, D. , 
and Cantor, c. , "Separation of Yeast Chromosome-sized" 
DNAs by Pulsed-field Gradient Gel Electrophoresis," 
cell 37:67, 1984.) By forcing the molecules to re-orient 
themselves between pulses, longer 'DNA fragments can be 
separated from shorter fragments. However, the induced 
electrical field is not uniform throughout the gel with 
this technique. The macromolecules tend to migrate 
through the gel with a mobility and trajectory that 
depend on where in the gel the samples are loaded. This 
position-dependent distortion, due to the non-uniformity 
of the electrical field, hinders lane-to-lane comparisons 
of the bands. 

Carle et al. describe a simple variant of PFGE, 
field-inversion gel electrophoresis (FIGE) , in which the' 
alternating fields are 180' apart. The polarity is 
reversed at intervals with a longer forward than reverse 
pulse time, thereby generating a net forward sample 
migration. This technique provides acceptable resolution 
up to 800 Kb (Kilobases) . For fragments above a certain 
size, however, large and small DNA molecules tend to 
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migrate together. By gradually changing (ramping) the 
voltage over the course of the electrophoresis, one can 
extend the separation range before this occurs. However, 
a set of parameters appropriate for the size of the 
molecule to separate must be selected at the outset. 
(See, Carle, G. , Frank, M. , and Olson, M. , 
"Electrophoretic Separations of Large DNA Molecules by 
Periodic Conversion of the Electric Field," Science 

232:65, 1986.) 

Gardiner and Patterson describe a transverse 
alternating field electrophoresis (TAFE) system in which 
the gel is oriented vertically and the usual four- 
electrode array is not placed in the plane of the gel but 
in front and in back of it. (Gardiner, K. , Laas, W., and 
15 • Patterson, D. , "Fractionation of Large Mammalian DNA 
Restriction Fragments Using Vertical Pulsed-Field 
Gradient Gel Electrophoresis," Som. Cell Molecular 
Genetics 12:185, 1986.) 

Chu et ai. teach that the critical factor in 
20 resolving large DNA fragments is keeping the angle - 

between the two directions of the alternating electrical 
field greater than 90*. (Chu, G. , Vollrath, D. , and 
Davis, R., "Separation of Large DNA Molecules by Contour- 
clamped Homogeneous Electric Fields," Science 234:1582, 
25 1986.) A contour-clamped homogeneous electric fields or 

CHEF system employs a hexagonal array with pairs of 
electrodes placed at 120'. Using transistors to force 
the electrodes to their ideal voltages, the field is 
clamped. Increased resolution is obtained without co- 
3 0 migration of samples. CHEF-DR II (BIO-RAD, Richmond, 

California) employs this technique and is used 
extensively for the separation of large DNA molecules. 
(U.S. Patent Number 4,878,008, Contour-clamped 
homogeneous electric field generator, issued to C. 
35 Ragsdale.) 

Current systems, however, require the user to 
manually determine and enter the parameters for the 
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electrophoretic operation. This technique is 
inefficient, requiring Much trial and error bv ^. 
specifically, this approach fails to t I 
pledge base that may be ^I^IT'' " ' 

:i :: p 0 : ret t c separation of dna . 

variety of set conditions. The invention ^ • , 
^.s an, apparatus which ^1;^™™ 

= P h=— — : ~ - 

aacromolecules with the aid of a computer Th* «. 
obtained by separating aacroaolecules in Dart , ° 

largest J" ° Per " i0n ' th * «".r. the smallest and 

largest dha sizes to be resolved. From this inform.n 
the system automatically calibr.ti.. , information, 

device with tk calibrates an electrophoresis 

switch 71 T, V ° lti9e ° radient ' An9le ' Init ^ 

witch Time, Final Switch Time. Forward Voltage Gradient 

and Reverse Voltage 0r .di.nt parameters, thereby 
increasing the resolution of electrophoresis. 

Ptasent J,^!'' " ^ ^ ° £ th ° ™ « «- 

th. el.ctro^:res B i: S d :v"r ^ " "«~ " 
Program. ' " * fl ° WChart ot th * ™EF H.pper 

Fig. 3 is a flowchart of case l. 

of the switch 3 ; "" 8 "* * tl0Webart ° £ the d «^ination 
oi -cne Switch Times for Case 1. 

Rcutin.. 5 " " "° WChart ° f thS B06111 ^ ^.raging 



Fig. 6 is a flowchart of Case 2. 

Fig. 7 is a flowchart of the determination of 
the initial switch Times for Case 2. 

Fig. 8 is a flowchart of the determination of 
the Final Switch Times for Case 2. 

Fig. 9 is a flowchart of the Mobility 
determination for Case 2 . 

Fig. 10 is a flowchart of Case 3. 

Figs. 11A-B are a flowchart of Case 4. 

Fig. 12 is a flowchart of the determination of 
the Initial Switch Times for Case 4. 



Overview 

The present invention provides methods and 
apparatus which work in conjunction with an 
electrophoresis means to improve the separation and 
purification of macromolecules, particularly DNA. 
According to the invention, a user enters information 
about the macromolecular sample to be separated. 
Parameters which optimize separation are then determined. 
Finally, the electrophoresis means is calibrated using 
these parameters. 

Referring to Fig. 1A, the present invention may 
be embodied in a system such as the system 100, which 
comprises a computer system 150, a CHEF Mapper 120, and 
an electrophoresis device or means 130. Computer system 
150 further comprises a central processor 108, a system 
memory 107, an input means or keyboard 103, a mouse 102, 
a fixed disk 104, a display monitor 101, an external 
interface 105, a printer 109, and an input/ output 
controller 106. The components of computer system 150 
are coupled to a system bus 115, thereby providing a link 
between all components. In a preferred embodiment, an 
appropriately programmed IBM compatible personal computer 
is used (International Business Machines, Armonk, NY) 
running under Microsoft windows (Microsoft, Redmond, WA) . 
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Referring to Fig. is, CHEF Mapper I20..and 
electrophoretic device 130 are further illustrated. 
Mapper 120 is coupled to device 130 and, thus can 
directly calibrate device 130 with optimized parameters. 
In a preferred embodiment, device 13 o is an 
electrophoresis means employing CHEF techniques. 
Macromolecule samples, which are initially placed in 
wells, migrate across a gel in response to an electric 
field. More specifically, the samples migrate under the 
influence of an electric field vector generated by a 
plurality of electrodes. 

In one embodiment, CHEF Mapper 12 0 is also 
coupled to external interface 105, such as a RS-232 
serial port, thereby allowing Mapper 120 to receive 
output from system 150 which, in turn, is used to 
calibrate electrophoresis device 130. m another 
embodiment, CHEF Mapper 120 is coupled to its own input 
device, such as a keyboard 121 and/or a bar wand (bar 
code reader) 122, thereby allowing Mapper 120 to receive 
information directly from the user, from alternative 
input formats, or from other systems. As shown in Fig. 
1A, Mapper 120 can also be coupled to external interface 
105 and to keyboard 121 and bar wand 122, all at the same 
time. 

System 100 operates in a fully automatic mode: 
the user enters the lowest and highest fragment sizes (in 
kilobases) of macromolecules (typically DNA) expected to 
be separated. A knowledge base of relationships 
discovered by applicants is provided which enables system 
100 to automatically determine the optimal parameters: 
Forward Voltage Gradient, Reverse Voltage Gradient, 
Voltage Gradient, Angle, Initial Switch Time, Final 
Switch Time, and Run Duration, as well as the best 
electrophoretic technique. Thus, maximum macromolecule 
35 separation can be achieved. 

In an alternative embodiment, system loo 
operates in a semi-automatic mode: the user can specify 
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additional parameters, including agarose type, agarose 
concentration, buffer type, buffer concentration, and 
temperature. These values are then used, in conjunction 
with the DNA sizes, to determine the optimal parameters. 

The individual parameters for system 100 will 
now be described in further detail. The Forward Voltage 
Gradient expresses the voltage differential per distance 
(V/cm) which is used during the forward pulse of FIGE. 
The Reverse Voltage Gradient, on the other hand, is the 
voltage differential per distance (V/cm) for the reverse 
pulse. The Voltage Gradient is simply the difference 
between the Forward and Reverse Voltage Gradients. 

Switch Time is the duration that 
electrophoretic device 130 remains at a particular state 
15 • (before switching again). Thus, this is the period of 
time that the system maintains a particular set of 
vectors (voltage gradient at a particular angle) . The 
Switch Time is ramped (varied) during the course of 
electrophoresis. Preferrably, the Switch Time is ramped 
20 in an increasing fashion. More preferrably, the Switch 

Time is ramped in a non-linear, increasing fashion. In a 
preferred embodiment, the Switch Time is ramped 
logarithmically. In an alternative embodiment, the 
Switch Time is ramped hyperbolically. In yet another 
25 embodiment, the user specifies the ramping by defining a 

curve on a coordinate system, for example, a cartesian 
coordinate system where X is the normalized Run Duration 
and Y is the normalized difference between Initial and 
Final Switch Times. Regardless of the ramping technique 
3 0 employed, it is preferable to divide the Switch Time into 

the Initial and the Final Switch Times indicating the 
upper and lower limits of the Switch Time ramp. 

Angle and Run Duration retain their usual 
connotations. The Angle is the angle maintained between 
3 5 the alternating fields and the sample. Historically, 

forward or down is defined as angle equal to 0*; for 
reverse or up, angle equals 180'. Finally, Run Duration 
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out by electrophoresis device 13 o. 

set of tJ lnC i 3 ^ ^ faaUiar Wlth 3 

set of terminology, three different program types are 

Zt " 1 Pr6ferred emb ° diment: 1) ««. 2) Two-state 
CHEF, and 3) Multi-state CHEF. This grouping aids the 
user, for example, in viewing run parameters. The FIG- 
has only f orvar d (angle - 0 -> and reverse (angle = la0 ") 
Pulses. m th. Multi-state CHEF or program type, forward 
occurs when angle eguals 0-, but there are all! of fills 
either positive or negative. For example, 60- to one ' 
side has a corresponding offset of -60- so that the net 
load is xn a forvard direction. Two-state CHEF is 

15 Inal ilar tl MUlti " State CHEF ™« it uses an included 

angle Thus, for so- and - 60 . angles, there would be an 
included angle of 120 «. 

relet' J 1 " 6 th " e Pr ° 9ram tyPeS haVe s ^^t/ S uper S et 
relatxonshxps. All of the program types can be described 
by the Multi-State CHEF program type. This program ty^ 
20 may include up to eight blocks, each blocJc ^ "« 

includxng up to fifteen states (vectors). The Two-state 
program type includes one block with two symetrical 
states. All values for each paramcter ^ ^ 

are xdentical except for angle. The angles are 
symetrical so that they sum to zero (forward direction). 

CHEF t. T E ^ 3lSO * SPeCial ° f the W^ti-state 

CHEF, it also consxsts of one block with two vectors. 

However, these vectors are not symetrical. One vector 
always has an angle of o- (forward), while the other has 
an angle of 180- (reverse). 

For standard conditions in the fully automatic 
mode, the following values are assumed for small DNA 
sxzes (< 2500 Kb) : 

35 

1) agarose type » "Low EEO" (industry 

standard) , 
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2) agarose concentration - 1.0%, 

3) buffer type - TBE (tris-borate EDTA) , 

4) buffer concentration - 0.5X TBE, and 

5) buffer temperature = 14 degrees Celsius. 

5 

The following standard conditions are assumed for larger 
DNA sizes) : 

1) agarose type - Low EEO, 

10 2) agarose concentration = 0.8%, 

3) buffer type = TAE ( tris-acetate EDTA), 

4) buffer concentration - 1.0X TAE , and 

5) buffer temperature = 14 degrees Celsius. 

15 , The knowledge base of system 100 is derived 

from calibration data obtained by separating DNA 
molecules of known size under a variety of planned 
conditions chosen to optimize separations. In a 
preferred embodiment, the knowledge base is embodied in 

20 constants, functions, and tables as part of a CHEF Mapper 

Program. The functions were obtained with polynomial and 
curvilinear regressions. The surfaces generated were 
fitted with a distance weighted moving least squares 
procedure performed on normalized data. (fifee, Lancaster, 

25 Peter, and Salkauskas, Curve and Surface Fitting: An 

Introduction, Academic Press, 1974; McLain, Drawing 
Contours from Arbitrary Data Points, Comput. J., 17:318- 
324, 1974.) Both the functions and the fitted surfaces 
are used for interpolation and in some cases for 

3 0 extrapolation. 

Effect of PWV Size 

The fully automatic operation of system 100 is 
divided into four different cases, the actual case being 
determined by the maximum and minimum DNA sizes entered 
35 by the user. In a preferred embodiment, three ranges of 

DNA sizes are defined: (1) 1-50 Kb, (2) 50-2200 Kb, and 
(3) greater than 2500 Kb. Each range has its own 
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particular set of macromolecular characteristic* which 
can be exploited to optimize separation. 

For DNA fragments of 1-50 Kb, for example, 
movement in a pulsed field is a logarithmic function of 
its size. This is in contrast to larger DNA fragments 
which have a mobility which is linear to their size. 
Fragments of different sizes do not respond identically 
to fields of different strengths. Larger fragments, 
«.g., 30-50 Kb, move relatively faster at 9 V/cm than do 
fragments of about 1-25 Kb. Different field strengths 
can be used to reduce the large spacing (distance between 
DNA bands) that occurs with the smaller bands (less than 
10 Kb). As a result, more uniform spacing occurs over 
the entire range. Thus, a preferred embodiment employs a 
field-inversion technique which uses different voltages 
for the forward and reverse fields for fragments below 50 
Kb to spread out fragments in a more uniform manner 
(compared to pulsed fields of identical voltages). 

DNA between 50 and 2500 Kb represent the range 
that was first described by Schwartz and Cantor. Most 
pulsed-f ield gel separations are -designed to separate DNA 
within this range. 

To separate large fragments of DNA, the field 
strength must be reduced to prevent the large molecules 
from getting tangled or broken in the gel. However, 
these low fields strengths necessitate extremely long run 
times, since the angle between the two fields influences 
the mobility of the DNA, run times may be decreased by 
reducing the angle between fields. 
30 System n p «, m i"P 

The operation of the system 100 will now be 
described in detail, rig. 2 is an overview of the CHEF 
Mapper Program 200, which is a set of procedures or steps 
stored in system memory 107 for instructing central 
35 processor 108. in another embodiment, CHEF Mapper 

Program 200 is stored in a firmware, e.g., a ROM, and 
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coupled either to system bus 115 or directly to Mapper 
120. 

To initiate program 200, in step 201 the user 
enters input parameters into keyboard 103 and/or mouse 
102; if program 2 00 is directly coupled to Mapper 120, 
parameters are more preferrably enter into keyboard 121 
or bar wand 122. For automatic mode, the user need only 
enter the minimum DNA size and the maximum DNA size of 
the sample to be separated. For the non-automatic mode, 
the user also enters the agarose type, agarose 
concentration, buffer type, buffer concentration, and 
temperature. Regardless of the mode, the input 
parameters are stored in a table, InputParms, located in 
system memory 107 and/or on fixed disk 104. 

In step 202, if the maximum DNA size entered 
("High") is less than or equal to 50 Kb, then the 
procedure goes to step 2 03 to perform Case 1. Otherwise, 
the procedure goes to step 204. If High is less than or 
equal to 2500 Kb in step 204, then the procedure goes to 
step 205 to perform Case 2. Otherwise, the procedure goes 
to step 206. If High is less than or equal to 6000 Kb in 
step 206, then the procedure goes to step 2 08 to 
determine if the minimum DNA size entered (Low) is less 
than or equal to 2000 Kb. If this is true, then the 
routine goes to perform Case 4 in step 210. otherwise, 
the procedure goes to perform Case 3 in step 209. 

However , if High is not less than or equal to 
6000 Kb in step 206, then in step 207 all output values 
are set equal to zero. After this step, or after Cases 1-4 
are performed, the procedure goes to step 211 and saves 
all the output values in an OutputParms table. In step 212, 
these values are sent to the CHEF Mapper 120, either 
directly by program 2 00 or via external interface 105. 
CHEF Mapper 120, in turn, calibrates electrophoresis 
device 130 with the optimum Angle, Initial Switch Time, 
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Final switch Time, Voltage, Run Duration, Forward 
Voltage, and Reverse Voltage. 
Individ!^] ?ffig1TS 

5 ?nn m w individual casM within CHEF Mapper Program 

200 which for* the knowledge base will now be described 
Fig. 3 illustrates the steps for Ca« i, that is, when High 
is less than or equal to 50 Kb. m step 301, field- 
inversion gel electrophoresis (FIGE) technique is chosen 

l\T* 3 ° 2 ' FOrWard V ° ltage G « di ««t is set equal to 

9 v/cm. in step 303, the Reverse Voltage Gradient is set 
equal to 6 V/cm. m step 304, the Angle is set equal to 
180- , and the Run Distance is set equal to 10 cm. In 
step 305, the switch Times are determined. 

FigS * 4A " B are a flowchart which illustrates 
the steps of determining the Switch Times, m step 401 
if Low is less than 25 Kb, then in step 402 the Initial ' 
Switch Time is set equal to 1.2 seconds. Otherwise, the 
procedure goes to step 403 to calculate the Initial 
Switch Time using the following relationship: 
20 (Function 1L) 

Y - 0.0 + 0.015797 X + 0.001365 X 2 

where X is the value of Low, and Y is the Switch Time in 
seconds. i„ step 404, if the Initial Switch Time is 
between 0.0 seconds and 0.05 seconds, then the Initial 
Switch Time i. .at equal to 0.05 seconds in step 405. 
Otherwise, step 405 is skipped. In step 406, the Final 
Switch Time i. determined from the following 
relationship: 
30 (Function 1U) 

Y « 0.0 + 0.009166 X + 0.000137 X 2 

where X is the value of Hi* and Y is the Switch Time in 
seconds. The result of Y i. multiplied by 1.15, an 
35 empirically derived adjustment factor, in step 407 if 

the Final Switch Tim. is between 0.0 seconds and 
0.05 seconds, then in step 408 the Final Switch Time is 
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set equal to 0.05 seconds. Otherwise, step 408 is 
sJcipped. In step 409, if the Initial Switch Time is 
greater than the Final Switch Time, then in step 410 the 
values for Initial and Final Switch Tines are reversed 
(swapped for one another). Otherwise, step 410 is 
skipped. The subroutine returns, thus concluding 
step 305. 

In step 306 if Low is between 1 Kb and 6 Kb or 
equal to either, then in step 308 the Mobility (cm/hr) is 
determined using the following relationship: 
(Function 2) 

Y » 0.65426181 - 0.0838514428 X + 0.0060879216 

X 2 , 



where X is the value for Low. However, if Low is not 
between the range of step 306, then in step 3 07 the 
Mobility (cm/hr) is determined from Table 1 (Appendix) 
and the Mobility Averaging Routine. 

Pig. 5 illustrates the steps of the Mobility 
Averaging Routine. In step 501, the Switch Time Range is 
determined from the Initial and Final Switch Times. In 
step 502, the number of iterations, W, is determined by 
dividing the Switch Time Range by five and rounding to 
the nearest integer. In step 503, the Delta Switch Time 
is determined by dividing the Switch Time Range by N. In 
step 504, the following constants are set: sum - 0.0, i 
- 0, Current Switch Time « Initial Switch Time. In 
step 505, i is tested to count N. If i is not less than 
or equal to N + l, then the routine divides sum by N in 
step 506 and stops, thus returning to step 3 07. 
Otherwise, the routine continues to step 507 to calculate 
a Mobility from either Table 1 or Table 2 (Appendix) 
using Low as the values of DNA size, current Switch Time 
as the Switch Time, and bilinear interpolation. In 
step 508, the Mobility is added to the sum. In step 509, 
the current Switch Time is incremented, in step 510, i 
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is incremented, and the procedure loops back to ..step 505. 

After completion of the Mobility Averaging 
Routine in step 307, the Run Time is determined in step 
309 by dividing the Run Distance by Mobility and 
multiplying by a Calibration Factor, slight variations 
in any one of several factors, e.g., water quality and 
laboratory technique, can significantly effect the 
operation of system 100. To minimize these variations, 
the user adjusts the Calibration Factor based on his or 
her experience with the system. The output of program 
200 is adjusted accordingly. The Cue i routine concludes 
by returning to step 203. 

Fig. 6 illustrates the steps of Case 2. in 
step 601, Angle is set equal to 120-, and Run Distance is 
set equal to 10 cm. In step 602, the Voltage Gradient is 
set equal to 6 V/cm. In step 603, Initial Switch Times 
for Case 2 are determined. Fig. 7 illustrates the steps 
involved, in step 701, if Low is less than 24 Kb, then in 
step 702 the Initial Switch Time is determined using the 
following relationship: 
(Function 22) 

Y - 0.0 + 0.072277 X - 0.008314 X 2 + 0.000581 X 3 



where X is Low, and Y is the Switch Time. The procedure 
then jumps to step 706. However, if Low is not less than 
24 Kb in step 701, then in step 703 Low is examined to see 
if it is less than 50 Kb. If this is true, then in 
«tep 704 the Initial Switch Time is determined using the 
following relationship: 
(Function 21U) 

Y - -0. 0701979 + 0.074324119 X + 0.00000733 6094 

X 2 
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where X is Low and Y is the Switch Time. The routine 
continues on to step 706, However, if Low is not less 
than 50 Kb in step 703, then in step 705 the Initial 
Switch Time is determined using the following 
relationship: 

(Function 21L) 

Y » 0.95604 + 0.11662974 X - 0.00000574 22 6 X*, 

where X is Low and Y is the Switch Time. In step 706, if 
the Initial Switch Time is between 0.0 seconds and 
0.05 seconds, then in step 707 the values for the Initial 
and Final Switch Times are inverted or reversed. The 
routine concludes, thus concluding step 603 of the 
calling routine. 

In step 604, the Final Switch Time is 
determined. Fig. 8 illustrates the steps of this method. 
In step 801, if High is less than 80 Kb, then in step 802 
the Final Switch Time is determined using the following 
relationship: 
(Function 21U) 

Y - -0.0701979 + 0.074324119 X + 0.000007336094 

X 2 

where X is High and Y is the Switch Time. However, if High 
is less than 80 Kb in step 801, then the procedure goes 
to step 803 to calculate the Final Switch Time using the 
following relationship: 

(Function 20) 

Y - 0.0 + 0.007605 X + 0.000686 X 2 , 

where X is High and Y is the Switch Time (sec) . The 
procedure continues on to step 804, where the Final 
Switch Time is multiplied by 1.15. In step 805, if the 
Final Switch Time is between 0.0 seconds and 
0.05 seconds, then in step 806 the values for the Initial 
and Final Switch Times are inverted. The procedure 
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concludes, thus concluding step 604 of the calling 
routine . 

In step 605, if the Initial Switch Time is 
greater than the Final Switch Time, then in step 606 the 
values for the Initial and Final Switch Times are 
inverted, otherwise, step 606 is skipped, m step 607, 
the Mobility for Case 2 is determined. 

Fig. 9 illustrates the process involved in step 
607. m step 901, if Low is less than 50 Kb, then in 
step 902 the Mobility (cm/hr) is determined using the 
following relationship: 
(Function 23) 

Y - 1.0124518 X • 0 *°"i222U 

where X is Low and Y is the Mobility, m step 903, if Low 
is less than 169 Kb and the Final Switch Time is less 
than lio.o, then in step 904 Mobility (cm/hr) is 
determined using Table 2 (Appendix) , Initial and Final 
Switch Times, Low as the DNA size, and the Mobility 
Averaging Routine (Fig. 5) . Then, the procedure 
concludes by returning. 

However, if step 903 is false, then the 
procedure goes to step 905 to test whether Low is less 
than 2500 Kb. If this is true, then in step 906 the 
Initial Switch Time is tested to see if it is between 
110.0 and 300.0 seconds. If this is true, then at 
step 907 the Mobility (cm/hr) is determined using Table 3 
(Appendix), the value of Low as the DNA size, and linear 
interpolation. Thereafter, the procedure concludes, thus 
concluding step 607 of the calling routine. 

However, if the known DNA size is not less than 
2500 Kb in step 905, then steps 906 and 907 are skipped. 
Furthermore, if step 906 is false, then step 907 is 
•kipped, in step 908, Low is tested to see whether it is 
l««s than 1169 Kb. If this is true, then in step 909 if 
the Initial Switch Time is less than lio.o seconds and 




the Final Switch Time is greater than 110.0 seconds, then 
the procedure goes to step 910 to calculate Mobility 
using Table 2 (Appendix) , Initial and Final Switch Times, 
Low as the DNA size, and the Mobility Averaging procedure. 
Then, the routine concludes by returning. 

However, if Low is not less than 1169 Kb in 
step 908, then steps 909 and 910 are skipped. Also, if 
step 909 is false, then step 910 is skipped. Under these 
circumstances, the routine arrives at step 911 where it 
sets all output values to zero before concluding 
step 607. 

In step 608, the Run Time is set equal to the 
Calibration Factor multiplied by the Run Distance and 
divided by the Mobility. The Case 2 procedure concludes 
and, thus, concludes step 205. 

Fig* 10 illustrates the steps for the Case 3 
procedure, which is called at step 209 (Fig. 2) . In 
step 1001, Angle is set equal to 106*. In step 1002, the 
Voltage Gradient is set equal to 2 V/cm. In step 1003, 
the Initial Switch Time is determined using the following 
equation: 
(Function 50) 

Y » 8.6 + 0.0053 X 

where X is High multiplied by 60 and Y is the Switch Time. 
In step 1004, if the Initial Switch Time is between 0.0 
and 0.05 seconds, then in step 1005, the Initial Switch 
Time is set equal to 0.05 seconds. Otherwise, step 1005 
is skipped. In step 1006, the Final Switch Time is set 
equal to the Initial Switch Time. In step 1007, the 
Mobility (cm/hr) is determined by using the following 
relationship: 
(Function 51) 

Y - 0.047 + 0.00000251 X + -1.56 10" 9 X 2 , 

where X is Low and Y is the Mobility. In step 1008, the 
Run Distance is set equal to 3 cm. In step 1009, the Run 
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Time is set equal to the Calibration Factor multiplied by 
the Run Distance and divided by Mobility. The Case 3 
concludes and, thus, concludes step 209. 

Figs. HA-B illustrate the steps of Case 4 
called at step 210 (Fig. 2). m step 1101, An gle is set 
equa to 106-. In step H02, the Voltage Gradient is set 
equal to 2 V/cm. m step 1103, the Run Time is set equal 
to the Calibration Factor multiplied by 72 hours, in 
step 1104, the initial Switch Time is set equal to 
20 minutes, m step 1105, the Final Switch Time is 
determined using the following relationship: 
(Function 50) 

Y = 8.6 + 0.0053 X 
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In 
15. 
0 and 



15 where x is High multiplied by 60 and Y is Mobility. 

step 1106, the Final Switch Time is multiplied by 1 15 
in step 1107, if the Final Switch Time is between 0.0 ana 

0. 05 seconds, then in step 1108, the Final Switch Time is 
set equal to 0.05 seconds. Otherwise, step 1108 is 
dipped, m step H09, the output values for block 1 are 
saved. 

in steps llio-iui, the Angle is set equal to 
120- and the Voltage Gradient is set equal to 6 V/cm. m 
step 1112, the Final Switch Time is set equal to 
3 minutes. In step 1113, Initial Switch Times are 
determined. 

Fig. 12 illustrates the method for determining 
the Initial Switch Times for Case 4. i„ step 1201, if Low 
is less than 24 Kb, then in step 1202 the initial Switch 
Time is determined using the following relationship: 
(Function 22) 

Y - 0.0 + 0.072277 X - 0.008314 X* + 0.000581 X 3 

9 

" where X is the value of Low and V is the Switch Time. The 

routine then continues on to step 1206. However, if Low 

1. not less than 24 Kb in step 1201, then the procedure 
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goes to step 1203 to test whether Low is less than 50 Kb. 
If this is true, then in step 1204 the Initial Switch 
Time is determined using the following relationship: 

(Function 21U) 

Y = -0.0701979 + 0.074324119 X + 0.00000733 6094 

X 3 , 

where X is the value of Low and Y is the Switch Time. The 
procedure continues on to step 1206. However, if Low is 
not less than 50 Kb in step 1203, then the Initial Switch 
Time is determined in step 1205 by using the following 
relationship: 

(Function 21L) 2 

Y - 0.95604 + 0.11662974 X - 0.00000574226 X , 



where X is the value of Low and Y is the Switch Time. 

The procedure continues on to step 12 06 to 
determine if the Initial Switch Time is between 0.0 and 
0.05 seconds. If this is true, then in step 1207 the 
Initial Switch Time is equal to 0.05 seconds. Otherwise, 
step 1207 is skipped. Step 1208 determines if the 
initial Switch Time is greater than the Final Switch 
Time. If this is true, then in step 1209 the values for 
the Initial and Final Switch Times are inverted. The 
25 routine concludes and, thus, concludes step 1113. 

in step 1114, the Mobility (cm/hr) is 
determined in the same manner as was described for Caie 2 
(Fig. 9). In step 1115, the Run Distance is set equal to 
4 cm. In step 1116, the Run Time is determined by 
30 multiplying the Calibration Factor by the Run Distance 

and dividing by the Mobility. The procedure concludes by 
returning. 

While the invention is described in some detail 
with specific reference to a single preferred embodiment 
35 and certain alternatives, there is no intent to limit the 
invention to that particular embodiment or those specific 
alternatives. The true scope of the invention is defined 
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not by the foregoing description but by the following 
claims. 
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Appendix 



Table 1. 



The first column (X) is the Minimum DNA Size [Low] 
(Kilobases) . 

The second column <Y) is the Switch Time (seconds) . 
10 The third column (Z) is the Electrophoretic Mobility 

(centimeters/hour) . 



Data 



15 . 


X 


Y 


Z 




8.14 


0.2 


0.273 




8.14 


0.5 


0.338 




8.14 


0.8 


0.384 


20 


10.18 


0.2 


0.245 




10.18 


0.5 


0.315 




10.18 


0.8 


0.362 




12.20 


0.2 


0.217 




12.20 


0.8 


0.351 


25 


15.00 


0.2 


0.180 




15.00 


0.5 


0.281 




15.00 


0.8 


0.348 




17.00 


0.2 


0.157 




17.00 


0.5 


0.262 


30 


17.00 


0.8 


0.339 




19.40 


0.2 


0.138 




19.40 


0.5 


0.236 




19.40 


0.8 


0.332 




22.60 


0.3 


0.177 


35 


22.60 


0.5 


0.216 




22.60 


0.8 


0.318 




24.80 


0.3 


0.158 
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24 . 80 


0.5 


0.191 


l4 • 80 


0 . 8 


0-306 


5Q on 


0.3 


0.128 


on 


0.5 


0.177 




0 . 8 


0.274 


j J . 50 


0.3 


0.144 


J J . 50 


0.5 


0.156 


33 . 50 


0.8 


0.245 


38.40 


0.3 


0.105 


38.40 


0.5 


0.126 


38.40 


0.8 


0.210 


48.50 


0.5 


0. 095 


48.50 


0.8 


0.146 
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The first column (X) is the Min 



imum DNA Size [Low] 



(Kilobases) . 

The second column (Y) is the Switch Time (seconds) 
The third column (Z) is the Electrophoretic Mobility. 



Fitted Surface 
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X 




Y 






Z 




8.1400e+000 


2. 


OOOOe- 


•001 


2 


•7129e- 


001 


8.1400e+000 


2. 


3000e- 


■001 


2 


•8031e- 


001 


8.1400e+000 


2. 


6000e- 


•001 


2 


-8722e- 


•001 


8.1400e+000 


2. 


9000e- 


•001 


2 


•9268e- 


001 


8.1400e+000 


3.2000e- 


•001 


2 


•9759e- 


001 


8. 1400e+000 


3. 


5000e- 


•001 


3 


-0269e- 


001 


8.1400e+000 


3. 


8000e-001 


3 


.0842e- 


001 


8.1400e+000 


4. 


lOOOe- 


•001 


3 


. 1486e- 


001 


8.1400e+000 


4. 


4000e-001 


3 


.2185e- 


001 


8.1400e+000 


4. 


7000e-001 


3 


•2911e- 


001 


8.1400e+000 


5. 


0000e-001 


3 


•3 611e- 


001 


8.1400e+000 


5. 


3000e- 


001 


3 


.4220e- 


001 


8.14006+000 


5. 


6000e- 


001 


3 


.4727©- 


001 


8.1400e+O00 


5. 


9000e- 


001 


3.5142e- 


001 



8.1400e+000 
8 .1400e+000 
8 .1400e+000 
8.1400e+000 
8.1400e+000 
8.1400e+000 
8.1400e+000 

1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.01586+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 
1.0158e+001 

1.2176e+001 
1.2176C+001 
1.21766+001 
1.2176e+001 
1.21766+001 
1.2176e+001 
1.2176e+001 



6.2000e-001 
6. 5000e-001 
6.8000e-001 
7.1000e-001 
7.4000e-001 
7.7000e-001 
8. 0000e-001 

2.0000e-001 
2.3000e-001 
2.6000e-001 
2.9000e-001 
3.2000e-001 
3.5000e-001 
3.8O00e-O01 
4.1000e-001 
4.4000e-001 
4.7000e-001 
5.0000e-001 
5.3000e-001 
5.6000e-001 
5.9000e-001 
6.2000e-001 
6.5000e-001 
6.8000e-001 
7.1000e-001 
7.4000e-001 
7.7000e-001 
8.0000e-001 

2.0000e-001 
2.30O0e-O01 
2.6000e-001 
2.9000e-001 
3.2000e-001 
3.5000e-001 
3.80006-001 



3 .5486e-001 
3 . 5794e-001 
3 . 6108e-001 
3.64*67e-001 
3.6903e-001 
3.7423e-001 
3.7958e-001 

2.4379e-001 
2. 5399e-001 
2. 6236e-001 
2. 6923e-001 
2.7539e-001 
2.8155e-001 
2.8814e-001 
2.95236-001 
3.0264e-001 
3. 1004e-001 
3.1720e-001 
3.2401e-001 
3. 3004e-001 
3.3519e-001 
3.3968e-001 
3.4383e-001 
3.4795e-001 
3.5227e-001 
3.5690e-001 
3.6177e-001 
3.6691e-001 

2.1715e-001 
2.2861e-001 
2.3852e-001 
2.4691e-001 
2.54356-001 
2.6154e-001 
2.6893e-001 
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1.2176e+001 
1.2176e+ 001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1.2176e+001 
1. 2176e+001 
1. 2176e+001 
1.2176e+001 
1. 2176e+001 



4* 1000e-001 
4. 4000e-001 
4 . 7000e-001 
5. 0000e-001 
5. 3000e-001 
5. 6000e-001 
5 rv 9000e-001 
6. 2000e-001 
6. 5000e-001 
6. 8000e-001 
7 . 1000e-001 
7 . 4000e-001 
7 . 7000e-001 
8 • 0000e-001 



2.7663e-001 
2.8453e-001 
2. 9237e-001 
2.9990e-001 
3 . 0694e-001 
3, 1339e-001 
3. 1930e-001 
3 . 2481e-001 
3. 3013e~001 
3, 3540e-001 
3 . 4071e-001 
3 . 4605e-001 
3 . 5132e-001 
3 . 5688e-001 



1. 4194e+001 
1.4194e+001 
1- 4194e+001 
1. 4194e+001 
l-4194e+001 
1. 4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 
1.4194e+001 



2 • 0000e-001 
2 . 3000e-001 
2 . 6000e-001 
2 . 9000e-001 
3 • 2000e-001 
3 .50QOe-001 
3 • 8000e-001 
4 • 1000e-001 
4 ,4000e-001 
4 .7000e-001 
5 •0000e-001 
5.3000e-001 
5.6000e-001 
5 -9000e-001 
6.2000e-001 
6 . 5000e-001 
6.8000e-001 
7.1000e-001 
7.4000e-G01 
7*7000e-001 



1. 9194e-001 
2.0479e-001 
2.1642e-001 
2.2635e-001 
2.3494e-001 
2 . 4292e-001 
2.5084e-001 
2.5893e-001 
2.6725e-001 
2.7564e-001 
2.8362e-001 
2.9064e-001 
2.9722e-001 
3 . 0374e-001 
3.1026e-001 
3.1678e-001 
3.2333e-001 
3 .2989e-001 
3.3643e-001 
3*4295e-001 
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1.4194e+001 


8. 


000Oe-OOl 


3 . 


4960e-001 


1.6212e+001 


2. 


0000e-OOl 


1. 


6937e-001 


1.6212e+001 


2. 


3000e-001 


1. 


8355e-001 


1.6212e+001 


2. 


6000e-001 


1. 


9689e-001 


1.6212e+001 


2. 


9000e-001 


2. 


0824e-001 


1. 6212e+001 


3 *2000e-001 


2 . 


1764e-001 


1.6212e+001 


3. 


5000e-001 


2. 


2594e-001 


1. 6212e+001 


3. 


8000e-001 


2 . 


3394e-001 


1.6212e+001 


4. 


1000e-001 


2 . 


4202e-001 


1.6212e+001 


4 . 


4000e-001 


2 . 


5029e-001 


1.6212e+001 


4 , 


7000e-001 


2. 


5870e-001 


1.6212e+001 


5 


. 0000e-001 


2 


,6677e-001 


1.6212e+001 


5 


. 3O00e-O01 


2 


.7406e-001 


1.6212e+001 


5 


. 6000e-001 


2 


.8112e-001 


1.6212e+001 


5 


. 9O00e-O01 


2 


.8842e-001 


1.6212e+001 


6 


.2000e-001 


2 


.9598e-001 


1. 6212e+O01 


6 


. 5O00e-O01 


3 


.0371e-001 


1.6212e+001 


6 


.8000e-ooi 


3 


.1154B-001 


1.6212e+001 


7 


. 1000e-001 


3 


.1938e-001 


1.6212e+001 


7 


. 4000e-O01 


3 


,2717e-001 


1.6212e+O01 


7 


. 7000e-001 


3 


.3486e-001 


1.6212e+O01 


8 


. OOOOe-001 


3 


.4240e-001 
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1.8230e+001 


2 . OOOOe- 


001 


1. 


5044e- 


001 


1.82306+001 


2.3000e- 


001 


1. 


.6570e- 


001 


1.8230e+001 


2 . 6000e- 


•001 


1 


. 8059e- 


•001 


1.8230e+001 


2 .9000e- 


■001 


1 


.9314e- 


•001 


1.8230e+001 


3.2000e- 


-001 


2 


.0283e- 


•001 


1.8230e+001 


3. 50006-001 


2 


. 1081e- 


-001 


1.8230e+001 


3.8000e- 


-001 


2 


.1834e-001 


1.8230e+001 


4 . lOOOe- 


-001 


2 


.2594e-001 


1.8230e+001 
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The first column (X) is the Minimum DNA Size [Low] 
(Kilobases) . 

The second column (Y) is the Switch Time (seconds). 
The third column (Z) is the Electrophoretic Mobility. 
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The first column (X) is the Minimum DNA Size [Lew] 
(Kilobases) . 

The second column (Y) is the Switch Time (seconds). 
The third column (Z) is the Electrophoretic Mobility. 
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Table 3 



column (X) is DNA Size (Kilobases) 
column (Y) is Mobility (cm/hr) 



X Y 

10 260 0.280 

335 0.269 

427 0.260 

581 0.243 

681 0.233 

15 774 0.219 

832 0.209 

913 0.200 

1051 0.173 

1600 0.129 

20 2200 0.100 

2500 0.085 



5 The first 

The second 
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CLAIMS 



1. A method for improving the seoaration of 
macromolecules in an electrophoresis system, with the aid 
of a computer comprising the steps of: 

entering a first size representing the smallest 
macromolecule to be resolved; 

entering a second size representing the largest 
macromolecule to be resolved; 

storing a knowledge base representative of 
calibration data obtained from separating macromolecules of 
known size under conditions which optimize separations; 

testing said first size and said second size against 
said knowledge base, thereby generating a set of optimized 
parameters ; and 

calibrating said system with said optimized 
parameters, thereby optimizing macromolecule separation. 

2. The method of claim l, wherein said testing step 

comprises: 

testing said first size and said second size against 
said knowledge base, thereby selecting an angle, a run 
distance, a voltage gradient, an initial switch time, a final 
switch time, and a run distance. 

3. The method of claim 1 or claim 2, wherein said 
calibrating step comprises: 

calibrating said system for optimal separation by 
using the angle, the run distance, the voltage gradient, the 
initial switch time, final switch time, and the run distance. 

4. The method of any one of claims 1, 2 or 3 
wherein said knowledge base storing step comprises: ' 



storing a first plurality of relations between 
macromolecule sizes and switch times; 
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storing a second plurality of relations between 
macromolecule sizes and mobilities; and 

storing a third plurality of relations between 
minimum macromolecule sizes, switch times and mobilities. 

5. An electrophoresis system for improving the 
separation of macromolecules comprising: 
an electrophoretic means; 

a computer having a memory and a processor, said 
computer being selectively coupled to said electrophoretic 

0eaI1S; input means, coupled to said computer, for entering 
a first size representing the smallest macromolecule to be 
resolved and entering a second size representing the largest 
macromolecule to be resolved; 

a knowledge base, coupled to said memory, for 
storing information representative of calibration data 
obtained from separating macromolecules of known size under 
conditions which optimize separations; 

programming means, coupled to said memory, for 
instructing said processor to test said first size and said 
second size against said knowledge base, thereby generating 
set of optimized parameters; and 

means for calibrating said electrophoresis means 
with said optimized parameters, thereby optimizing 
macromo 1 ecu le separ at i on . 

6. The system of claim 5, wherein said programming 

means comprises: 

programming means, coupled to said memory, for 
instructing said processor to test said first size and said 
second size against said knowledge base, thereby 
angle, a run distance, a voltage gradient, an initial switch 
time, a final switch time, and a run distance. 

7. The system of claims or claim 6, wherein said means 
for calibrating comprises: 



60 



means for calibrating said system for optimal 
separation by using the angle, the run distance, the voltage 
gradient, the initial switch time, final switch time, and the 
run distance. 

8. The system of any one of claims 5, 6 or 7 
wherein said knowledge base comprises: 

storage means, coupled to said memory, for storing a 
first plurality of relations between macromolecule sizes and 
switch times, a second plurality of relations between 
macromolecule sizes and mobilities, and a third plurality of 
relations between minimum macromolecule sizes, switch times 
and mobilities. 

9. The system of any one of claims 5, 6, 7 or -8 
further comprising Mapper means for receiving said optimised 
parameters, said Mapper means being selectively coupled to 
said computer through an external interface. 

10. The system of any one of claims 5, 6, 7 8 or 
device 10 Said eleCtr ° Ph0retic i- a CHEF electrophoretic 
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